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also discuss that the singularity of the Spinless Salpeter wavefunction at the origin is 
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1 Introduction 



A modern method to calculate a non-relativistic bound state problem is based on the Effective 
Field Theory (NRQED) which was proposed by Caswell and Lapagefl]. This is the advantage 
of this method that provides a set of systematic rules (power-counting) that allows an easy 
identification of all terms that contribute to a certain order in the bound state calculations. 

In NRQED calculation, relativistic kinetic vertex needs more attention respect to the 
other interactions [2]. In this paper, we consider this vertex in the unperturbed part of the 
Hamiltonian as well as the Coulomb interaction to obtain the positronium energy correction. 
For this purpose we derive the Spinless Salpeter (SS) equation for the positronium with Coulomb 
potential using the NRQED Lagrangian in Sec. 2. The SS eigenfunctions and eigenvalues in 
Sec. 3 and 4 are used to describe the contribution of the relativistic vertex correction to the 
energy spectrum of the positronium at the order of a 4 and a 6 . 



There are many equations incorporating relativistic kinematics which all have the same non- 
relativistic limit [3]. However, if an expansion is made in powers of the momentum, the various 
equations differ in terms next to leading order. Furthermore, most of them are phenomenolog- 
ical equations and which are not derivable from field theory without making gross approxima- 
tions. An approach which seems promising is the Spinless Salpeter equation to obtain spectrum 
of a two-body system. This equation, with appropriate assumptions, can be derived from the 
Bethe-Salpeter equation[4] which is manifestly covariant and can be obtained directly from 
quantum field theory. We now consider the Spinless Salpeter Equation [5] as follows: 



2 Spinless Salpeter Equation 




(i) 



where for positronium, in the center of mass frame, one has 




(2) 



or 




(3) 
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where m is the mass of the particle and V(r) is the time component of a 4- vector potential where 
in the Coulomb case V(r) = I y L . In the non-relativistic limit, Eq.(3) leads to a Schrodinger 
equation with the effective potential V. In fact this equation gives all relativistic kinetic correc- 
tions to the Schrodinger equation with such effective potential. However we can derive Eq.(3) 
using NRQED Lagrangian. The momentum space equation of motion for an off-shell, time 
independent e + e - e + e~ four point function in the center-of-mass frame is: 



E 



G(p,q;s) + J ^V(p,q')G(q,q';s) = (2 tt) 3 6®{p - q) , 



where 



(2 7T) 



E = y/s — 2 m e , 



(4) 



(5) 



is the center-of-mass energy relative to the electron-positron threshold and V(p, q), the poten- 
tial, is defined as 



V(p,q) = V C oui(p,q)+V BF (p,q) + V 4 (p,q) + V 4dei (p,q)+SH kiI1 (p,q), 



(6) 



introduced in Ref[2]. Considering 5H kin (p) to all orders will lead to the following equation: 

G(p, q, s) + (V Coal + V')G(p, q, s) = (2 tt) 3 5® (p - q) , (7) 
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or 



where 



2 x /p^+m} - 1 G(p, q; s) + (V Coul + V')G(p, q; s) = (2 tt) 3 5^ (p-q), 



(8) 



(V Coul + V')G(p, q;s)= I p(- V Coul (p, q) G(q, q', s) + [ ^ V'(p, q') G(q, q', s) , (9) 



(2 tt)* 



(2 tt)* 



and 



V'(p,q) = V BF (p, 0) + V*{p, q) + V* A «(p, q) ■ 



(10) 



Now equation (8) can be rewritten as, 



2 J p 2 + m 2_ { ^- s _ Vconi) G(p,q;s) + V'G(p,q;s)= (2 5^\p - q) . (11) 



Since the Coulomb potential and the relativistic kinetic vertices are the only interactions that 
must be treated exactly, we can ignore V in the above equation which leads to the Eq.(3) 
in momentum space. Therefore one can start with the wavefunction of the SS equation with 
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Coulomb potential and determines the corrections by perturbation as is usual in NRQED. 
Whereas the effect of relativistic kinetic vertices is in the SS wavefunction of Eq.(3). To this 
end by substitution \I/ = ^ in the radial part of Eq.(3) one obtains: 

+ I _ (l + «)(* + « + l) 1\ u( ) _ Q (12) 

where 



k = -(£ + 1/2) + yJ(e+l/2) 2 -a 2 /4 & p = Ear/A. 
The eigenvalues and eigenfunctions of Eq.(12) are respectively [6] [7]: 

„2 n-1/2 



E nJ = 2m 



i+ " 



4(n + k) 2 



(13) 



UZ+ K = N(—) T-H2u + 2k + 1) p^exp(-^-), (14) 
in which TV is an appropriate normalization constant. We can now expand -J^ 2 - in terms of a 

^ Sch 

(when n — 1, £ — 0) as 

^ss(r) = * Sc fc(r){l + a 2 ( - ^ 7 + - i In (mm)) +...}, (15) 

where 7 is the Euler number and *&sch( r ) is the ground state wavefunction for the Schrodinger 
equation with Coulomb potential, 



1 /ma\ 3 / 2 
^Sch{r) = -j= y— J exp (-mar/2). 



(16) 



The Fourier transformation of Eq.(15) is 

^ss(p) = Hp)sch + «V(P) + • • • , (17) 

where V\Sch(p) an d V?(p) are 

, , , 87r 1 / 2 (ma/2) 5 / 2 , x 

*>&>> = - r/'^T^ {arctan(^)^_ (ma/2 )A 
{(ma/2) z + p z yp {. ma \ J 

+p(ma/2) In f ( ma / 2 ) 2 +P 2 \ + 3p ( mQ; /2)}. (19) 

We are now ready to calculate positronium energy corrections at the order of a 4 and a 6 , using 
SS wavefunction, Eq.(15). 
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3 a 4 Positronium Energy Corrections 

The known vertices in NRQED Lagrangian are given in [11]. The power-counting of NRQED[12] 
specifies all the diagrams contributing to the energy shift at the order of a 4 . These diagrams 
are all shown in Fig.(l) which can be calculated [8] to give the well-known positronium energy 
correction at the order of a 4 [9]. Contribution of relativistic kinetic corrections at this order are 
shown in Figs. (l(i,j)) which reads: 
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Figure 1: The whole bound state diagrams contributing to positronium energy corrections at 
the order of a 4 , ip is schrodinger wavefunction and C4 is the leading order result coming from 
tree-level matching (c 4 = cf^ = —™). 
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16n 3 (£+l/2) 64n 4 

To obtain positronium energy correction at the order of a 4 using SS equation, we should omit 
Figs.(l(i,j)) and calculate the remaining diagrams in Fig.(l) using SS wavefunction. That is 
obviously leads to the previous result at the order of a 4 and also extra pieces which start at 
higher order. These terms are irrelevant to the order of our interest. Meanwhile the result of 
Eq.(20) can be easily obtained by expanding the Eq.(13) in terms of a, 

1 _1 / 2 _ ma 2 ma 4 3ma 4 



E„ e = 2m 



1+ ° 



2m 4n 2 16n 3 (£+l/2) + 64n 4 + "" (21) 



4(n + k) 

4 a 6 Positronium Energy Correction 

The full calculation of the positronium energy correction at the order of a 6 , using NRQED has 
been done in reference [10]. In this paper, we focus on the positronium Hyperfine Splitting 
(HFS) at the order of a 6 coming from one-photon annihilation [14] [2]. To be more specific in 
the way that the SS wavefunction can be useful in the bound state calculation, we first briefly 
review the NRQED method in the following subsection. 



4.1 Matching and bound state energy shift at NNLO 

Since the single photon annihilation of electron-positron occurs in S — 1 state, one should 
calculate the diagrams which contain spin-1 Four-Fermion Vertex, 04(1—7 ann). We can write 
04(1 — 7 ann) = cf\l — 7 ann) + c^\l — 7 ann) + . . . where the superscript " (0)" and " (1)" 
indicate that the coefficients c 4 (l — 7 ann) have been derived using tree-level and one- loop 
matching, respectively. To obtain the contribution of HFS at the order a 6 , we should perform 

(2) 

matching at two- loop level to get c\ (1 — 7 ann). For this purpose, it is more convenient to 
write: 

c\ '(I - 7 ann) = c\ '(I - 7 ann) PC + c\ '(I - 7 ann) V c, (22) 

where PC stands for contribution coming from propagator corrections while VC for vertex 
corrections. Each terms in right hand-side of Eq.(22) can be obtained from Figs. (2-3). All 
the NRQED diagrams which contribute to HFS at the order of a 6 can be identified using 
the NRQED power-counting rules[12]. These are shown in Figs. (4-5) which are completely 
calculated in [8] [2]. Diagrams of Figs.(4(k,l,m)) which contain relativistic vertex corrections 
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Figure 2: The matching which determines c\' (1 — 7 ann) PC while c 4 (1 — 7 ann) PC can be 
similarly obtained from the one-loop matching procedure. 

result in [8]: 

. ^ 1 ,_,,„,. „ , . \ ma 6 , / A \ 1 fi ma 5 A . „. 

Fig{A, k) + Fi</ 4, 1) + Fi<7 4, m = — - In + -ma 6 + . (23) 

8 \maj 62 An m 

4.2 NNLO bound state energy shift using SS wavefunction 

As shown in section 2, the whole relativistic kinetic vertex as well as Coulomb vertex can be 
considered in unperturbed part of the Hamiltonian, Eq.(ll). Therefore, one should expect to 
omit relativistic kinetic vertex and using SS wavefunction instead of Schrodinger wavefunction, 
in all order of perturbation, to obtain the same results within the framework of NRQED. That 
is basically reduce the number of counter terms involve in this calculation [15]. Since our goal 
is to obtain the final result using SS wavefunction, we should omit Figs.(4(k,l,m)) and calculate 
the remaining diagrams in Fig(4-5) with the SS wavefunction. Straightforward calculations 
shows that we have the previous result for these diagrams and also extra pieces at higher order 
than a 6 . At first glance, it seems that there is no way to get the value of Figs.(4k,l,m) which 
we have omitted, but if we consider the diagrams which contributed to HFS at the order of a 4 
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Figure 3: The matching which determines c 4 (1 — 7 ann)vc while c 4 (1 — 7 ann) vc can be 
similarly obtained from the one-loop matching procedure. 



7 



Four-Fermi ~ T ann ) + c i 2) (! - 7 ann) ^ V>C ^ 



Spin-Spin ^ 



Transverse^ 



Darwin 



Relativistic , 
Kinetic ^ 



Derivative 




(k) 



1p tjj 



1p 1p 



VC / \ 
(a) 




^ -0 




Figure 4: All the bound state diagrams which contribute to a 6 except the diagrams with 
the Double Annihilation Interaction, is called the Derivative interaction comes from the 
Taylor expansion of each vertex of the one photon annihilation diagram as well as the photon 
propagator correction which is equal to = ^r- 

we get some pieces at the higher order. That is basically due to the fact that we should replace 
ipSch with ipss i n Fig-(l>k) and also using the p/m expansion of relativistic propagator. This 
means that the Fig.(l,k) should be replaced by the left hand side of Fig. (6). It is easy to show 
that Fig. (6(a)) leads to the previous result at the order a 4 , while the Figs.(6(c,e,f)) contributes 
to the higher order than a 6 . The only remaining diagrams relevant to our calculation are the 
Figs.(6(b,d)): 
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Fig(6(d)) = 
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(c) (d) 

Figure 5: The whole Double Annihilation bound state diagrams. 

The sum of Fig. (6(b)) and Fig.(6(d)) is just equal to (23). It should be noted here that 
the divergence in the first term of (24) is a direct consequence of the singularity of the SS 
wavefunction (15). In fact in the Schrodinger case Figs. (2(1, m)) are responsible to cancel the 
logarithmic divergence coming from the other diagrams while in the SS case this is due to the 
SS wavefunction. 




Summary 

In this paper, we have shown that the spinless Salpeter equation (3) correctly predicts the 
energy spectrum of a non-relativistic two body system such as positronium. In this way we 
have ignored the relativistic vertex correction in the bound state and therefore the calculations 
are made easier. On the other hand, although the SS wavefunction is singular at the origin (see 
Eq.(15)) but this is a crucial need to cancel UV divergence in the bound state system, Eq.(24). 



Acknowledgement 

The authors gratefully acknowledge the financial support of Shiraz University research council 
and IPM. 



9 




Figure 6: The diagrams at the order a 6 and higher coming from one photon annihilation at the 
order a 4 . 
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